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Abstract 

The production rates and transverse momentum distributions of the strange- 
beauty mesons B s and B* at pp colliders are calculated assuming fragmenta¬ 
tion is the dominant process. Results are given for the Tevatron in the large 
transverse momentum region, where fragmentation is expected to be most 
important. 


I. INTRODUCTION 

Data on the production of strange B mesons is rapidly accumulating. At LEP the DEL¬ 
PHI Collaboration has measured the probability for a b antiquark to hadronize into a weakly 
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decaying strange B meson |L|. Assuming the hadronization is dominated by the fragmen¬ 
tation of the heavy b antiquark, we determined the fragmentation functions Di_> Bs (z, p) 
and D^ B *(z, p) from the measured total probability for a b antiquark to hadronize into 
the lowest strange-beauty states B^So) and B*( 3 Si) ]2|. Here 2 is the fraction of the b 
antiquark momentum carried by the B s or B* at the scale p. The momentum distributions 
of the B s and B* mesons produced in e + e~ annihilation at the energy of the Z mass were 
then predicted |2|]. 

Here we extend our calculations to the production of the B s and B* states in pp an¬ 
nihilation at the Tevatron, using the fragmentation functions previously determined from 
the e + e _ data taken at LEP. Sec. II presents the cross section for the process in terms of 
the parton distribution functions and subprocesses and the fragmentation functions, which 
we discuss in Secs. Ill and IV, respectively. In Sec. V we present numerical results for 
the transverse momentum (px) distributions of the strange B s and B* mesons, as well as 
the total cross section for px > p™ in . Finally, we discuss these results and summarize our 
conclusions. 


II. CROSS SECTIONS 


In the parton model the cross section for the production by fragmentation of strange B 
mesons in proton-antiproton annihilation involves three main factors: the structure functions 
of the initial partons in the proton and antiproton, the subprocess in which the initial partons 
produce a particular parton in the final state, and the fragmentation of this final parton into 
the strange B meson. The transverse momentum distribution of the B s meson is then of the 


form 

da 

dp T 



Here fi/ p (xi ) and fj/p{xj ) are the structure functions of the initial partons i and j carrying 
fractions of the total momenta Xi and Xj in the proton and antiproton, respectively. The 
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production of parton k with momentum pt/z is described by the hard subprocess cross 
section & at the scale p. And D k ^ Ba (z: p) is the fragmentation function for the parton k 
to yield the meson B s with momentum fraction z at the scale p. The factorization scale 
p is of the order of the transverse momentum of the fragmenting parton. The soft physics 
is contained in the structure functions, while the parton subprocess is a hard process and 
can be reliably calculated in perturbation theory. We shall also assume the fragmentation 
of b into B s is also a hard process, an assumption that will require further discussion below, 
since the s quark is not very heavy. 

Clearly, there is a formula very similar to Eq. ([!]) for the production of the B*, which 
we need not write explicitly. Because the radiative decay B* —> B s + 7 is so fast and the 
B* — B s mass difference is so small, the transverse momentum of B* is hardly affected in 
the decay B* —> B s + 7 , and thus contributes to the inclusive production of the B s . We 
will, therefore, present each individually in the figures showing da/dpx and u{jpt > p™ in ) 
and also include their sum in a table, for convenience. 

In the fragmentation process k —> B s we have included gluons as well as b antiquarks. 
Although the direct g —> B s fragmentation process does not occur until order there 
is a significant contribution coming from the Altarelli-Parisi evolution of the b antiquark 
fragmentation functions from the heavy quark mass to the collider energy scale Q, which is 
of order a^ log (Q/m,b) due to the splitting g 6 . This has been shown to be significant for 
the production of B c mesons [[| and we have also included it here. 

III. STRUCTURE FUNCTIONS AND PARTON SUBPROCESSES 

In Eq. ([[]) we have assumed a factorization scale /i, with the parton distribution func¬ 
tions being due to soft processes below this scale, while the parton cross sections for the 
subprocesses are due to hard contributions above this scale. For the parton distribution 
functions we used the most recent CTEQ version 3 ||, in which we chose the leading-order 
fit, since our calculation is also a leading-order one. For the subprocesses we used the tree- 
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level cross sections, to be consistent, since the parton fragmentation functions, which are 
discussed in the following section, were calculated only to leading order. For the production 
of b antiquarks we included the processes gg bb, gb —> gb, and qq —> bb, while for the 
production of gluons g we included the processes gg —► gg, gq(q) —■> gq(q), and qq —> gg, 
where q denotes any of the quarks, u, d, s, c, b. 

For the running strong coupling constant a s we used the simple one-loop result evolved 
from its value sX g = m Z - 


a s (g) 


_ a s (mz) _ 

1 + ^f as (m z ) log (£) 


( 2 ) 


Here rif is the number of active flavors at the scale g and we chose a s (m z ) = 0.117 0. From 
the experimental results on 6-quark production cross sections at the Tevatron we know that 
there are discrepancies between the calculations and the experiments about a factor of two. 
Therefore, adopting a more complicated form for the strong coupling constant, or using the 
next-to-leading order subprocess cross sections, or taking the next-to-leading order structure 
functions is not critical. At this point, we are primarily interested in how important the 
fragmentation process is for B s and B* production and how the shapes of the calculated 
momentum distributions compare with forthcoming data. 

The factorization scale g in Eq. ([!]) deserves further discussion before proceeding. The 
scale g superficially appears to have been invented only to separate the production of B s 
mesons into structure functions, subprocess cross sections, and fragmentation functions. In 
principle, this is so and the production is independent of the choice of g. But in practice, 
the results do depend on g, since only if all the factors are calculated to all orders in a s 
will the dependence on g cancel. However, for example, the fragmentation functions were 
calculated only to leading order and, consequently, the results will depend on g. We shall 
investigate the sensitivity of the results to the choice of g by varying g from g R j2 to 2 g R , 
where g R = (parton) + ml is our central choice of the scale g. 


4 





IV. FRAGMENTATION FUNCTIONS 


To obtain the fragmentation functions at the scale y we have numerically integrated the 
Altarelli-Parisi evolution equations || from the scale y o, which is of the order of the 6-quark 
mass. The evolution equations for the fragmentation functions D b ^ Bs (z, y) and D g ^ Ba (z, y) 
are 


0 /* 1 dy r 1 dy 

y~ B r D b^ Bs (z,fi)= —Pb^b(z/y,v)D- b ^ B (y,n)+ —P b (z/y 1 y)D g ^ Ba (y 1 y), (3) 

oy Jz y J z y 

and 

d /* i dy f ^ dy 

y^-D g ^ Bs (z,y)= —Pg^- b (z/y,y)D- b ^ Bs (y,y)+ —Pg^g{ z /y,y)Dg^ Bs (y,y), (4) 

cj i~l j z y J z y 

with similar, coupled equations for the B* fragmentation functions. The splitting functions 
Pi^j{x, y), to leading order in a s (y), are explicitly given in Ref. 0. As boundary conditions 
in these calculations we have used the fragmentation functions at the scale /x 0 calculated in 
Ref. 0: 


Pb^B s ( z i yo) 


2a s (2m s )-\R(0)\- rz(l — z) 2 
8 l7rm^ (1 — (1 — r)z) 6 
x [6 - 18(1 - 2r)z + (21 - 74r + 68 r 2 )z 2 


-2(1 - r)(6 - 19r + 18r 2 )^ 3 + 3(1 - r) 2 (l - 2r + 2r 2 )^ 4 ] . 


( 5 ) 


Here z is the fraction of the b antiquark momentum carried by the B s meson, r = m s /[mb + 
m s ), and i?(0) is the B s meson A-wave radial wavefunction at the origin. The corresponding 
fragmentation function for b —> B* is [7] 


Dr 


b—>BZ ( z i ^0) ~ 


2 ay(2m s ) 2 |i?(0)| 2 rz(l — z)^ 


27-Krnl 


(1 — (1 — r)z) 6 


x [2 - 2(3 - 2 r)z + 3(3 - 2r + 4r 2 )z 2 
—2(1 — r)(4 — r + 2 r 2 )z 3 + (1 — r) 2 (3 — 2r + 2r 2 )z 4 ] , 


( 6 ) 


The induced gluon fragmentation functions D g ^ Ba)B *(z, y) are of order c^log^y/m b ) and 
become important at large values of the scale y relative to the b antiquark fragmentation 
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functions Di^ Bs B *(z,/d), which are of order a 2 s (n). The boundary conditions for the gluon 
fragmentation functions are 


D g ^B s {z, aO — D g ^B*{z, +) — 0 


( 7 ) 


for < 2(mi, + m s ), the threshold for producing B s or B* mesons from a gluon. 

The value of |i?(0)| J can be determined from the decay constant Jb s of the B s meson 
using the relation || 


2 = 3\mi 

lBs vr M Bs 

with f Bs = 207 ± 34 ± 22 MeV 0 and M Bs = 5.375 GeV @. 

If we choose the 6 quark mass to be mb = 5 GeV the only remaining parameter in these 
fragmentation functions is m s . This strange-quark mass parameter rn s has been determined 
previously 0 using the same initial fragmentation functions from the experimental value 
0 of the probability ff for a b antiquark to fragment into weakly decaying strange-beauty 
mesons: 


f: = / dz 

Jo 


Db^B a ( z > +o) + D b ^ B * (z, /+) 


( 9 ) 


Since the total probability for the b antiquark to fragment into a B meson is independent 

of scale, the initial scale /j,q in Eq. 0), which is of the order of the 6-quark mass, was chosen 

+47 

to be fio = mb + 2 m s as in Ref. [|7|. We previously found |m s = 318—24 MeV/c 2 based on 
the DELPHI measurement ff = 0.19 ± 0.06 ± 0.08 and assuming a s (mz ) = 0.12. Prior to 


the DELPHI measurement there were older data from LEP 10 which were consistent with 


this measurement. There is now also a new measurement from CDF |H[ of the quantity 

a(B s ) 


cr (B u ) + a(B d ) 


+.17 

= 0.26—.08 ± .08 


( 10 ) 


Since the total fragmentation probability P(b —» baryon) « 5 — 10%, the corresponding value 
of /“ = P(b — > B SJ B *) can be estimated assuming P(b —> B u +B d +B s +ba,rjons) = 1. Using 
the newer CDF data, Eq. (IToj), we obtain /“ = 0.196(0.186) for P(b —> baryon) = 5(10)%, 
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which is consistent with the DELPHI measurement. In the present calculations, therefore, we 
used the same value = 0.19±0.06±0.08 JIJ but a slightly different value of a s (m z ) = 0.117 

|J. For the strange quark mass parameter we then obtain 

+47 

m s = 298-23 MeV (11) 

using the scale-independent relation Eq. (||) . This value of m s is not significantly different 
from the previous determination and we shall use it in our fragmentation functions Eqs. (0) 
- ([|). The strong coupling constant in Eqs. (H) - © then has the value a s (2m s ) = 0.745. 
While this value of a s is uncomfortably large we note that it only enters, together with m s 
and -R(O), in the normalization of the fragmentation functions, which has been determined 
empirically through Eq. (|9j). 

V. RESULTS AND DISCUSSIONS 

The evolution equations were numerically integrated and the fragmentation functions at 
the scale Pr = \Jpj, (parton) + m 2 were combined with the structure functions and parton 
cross sections for the subprocesses to obtain the cross section, Eq. ©■ Specifically, we 
calculated the transverse momentum distributions of B s and B* mesons produced in pp 
collisions at the Tevatron energy y/s = 1.8 TeV. We assumed a cut-off on the transverse 
momentum of px > 5 GeV/c 2 and considered only the rapidity range \y\ < 1. In Fig. ]l] we 
show the cross section da/dpx for both B s and B*. To investigate the sensitivity of these 
results to the scale p we have also included the results for p = 2 p R and p = Pr/2. When the 
scale p is less than p 0 = rrib + 2 m s , which only happens for the case of p = pr/2 , we chose 
the larger of (p,po). From Fig. [I] it is clear that the choice of scale p is not critical for the 
transverse momentum distribution; in fact, the variation over the range Pr/2 < p < 2 pr is 
comparable to the current discrepancies between the measured and calculated b production 
cross sections. As one might expect, the ( 3 S'i) B* cross section is larger than the ( 1 S'o) B s 
cross section at all pr, however, their ratio is not precisely the naive prediction 3:1, but is 
about 20% smaller. 
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In Fig. we show the total cross section for the production of B s and B* mesons with 
transverse momentum above a minimum value As in Fig. [I] only the range p B > 5 GeV 
and \y\ < 1 was considered. The sensitivity to the scale p was investigated, as before, by 
considering p = Pr/2 and p = 2pR, and the choice of this scale is clearly not critical. Figure 
provides useful estimates of the production rates of B s and B* mesons at the Tevatron due 
to the fragmentation process. For convenience, we have also presented these cross section 
estimates in Table p} including the sum a(B s ) + cr(B*). 

The production rates and transverse momentum distributions for B s and B* meson pro¬ 
duction in pp collisions presented here assume that fragmentation is the dominant process 
and, therefore, are probably underestimates. At large enough transverse momentum the 
fragmentation process should dominate, as it falls off more slowly, even though it is only a 
part of the full order aj contribution. We also note that in our calculation the light quark 
fragmentation contribution is not included, since it is of even higher order in a s than the 
induced gluon fragmentation. At the initial heavy quark mass scale, the light quark frag¬ 
mentation function D q ^ Bs (z, p) is of order a A s (af) for q = u,d(q = s). But the D s ^ Bs (z, p) 
has been shown [7] to be suppressed by {m s /m b ) 3 relative to D b ^ Bs (z, p). Therefore, at the 
initial scale none of these light quark fragmentation functions are important. The light quark 
fragmentation functions might acquire a logarithmic enhancement from the evolution of the 
Altarelli-Parisi equations, but the splitting kernels P g ^ q (x, p) and P q ^ g (x, p) are of order a s , 
while P b ^ q (x,p) and P q ^b{x,p) are of order a^. Therefore, the most important sources of 
induced light quark fragmentation come from (i) P q ^i0 D b ^ Bs (z) , and (ii) P q ^ g ® D g ^ Bs {z). 
Hence, the order of the induced light-quark fragmentation functions at the scale p can be 
at most of order a f x log(p/m b ), which is suppressed by a s relative to the induced gluon 
fragmentation. Since the induced gluon fragmentation only contributes at about the 15-20% 
level, which is mostly due to the large gluon luminosity at low x, the induced light-quark 
fragmentation contribution is expected to be very small, at the level of a few %. And thus, 
the light quark contribution can be safely ignored in compare with other uncertainties in 
the calculation. 
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Our fragmentation functions, which we used for the boundary conditions for the Altarelli- 
Parisi equations at the scale of the heavy b quark mass, are rigorously correct in perturbative 
QCD if the constituent quarks are much heavier than Aqcd; e.g, B c mesons. Also higher 
order QCD corrections can, in principle, be systematically calculated. In the case of the B s 
meson, higher order corrections in the fragmentation functions might be large due to the 
fact that the strong coupling constant is evaluated at 2 m s and is equal to about 0.75. But 
we fitted the parameter m s to the experimental measurement of the total probability for 
b —> B s + B* and m s turned out to be about 300 MeV/c 2 . We can, therefore, be confident 
of the overall normalization, which is proportional to a^(2m s )\R(0)\ 2 /m\. The shapes of the 
fragmentation functions are only moderately sensitive to the value of rn s and comparing our 
results for the momentum distribution will, therefore, primarily test the assumption that 
their shapes are adequately given by perturbative QCD at the scale of the b quark mass. 

Our approach also appears somewhat similar to the approach used in some Monte Carlo 
programs, e.g., PYTHIA |TI| and HERWIG [13|]. However, the main difference lies in the 
fragmentation of the partons. The fragmentation models used in PYTHIA or HERWIG 
are rather different in spirit from our perturbative QCD fragmentation, and is based on 
string fragmentation, though the actual fragmentation models used in these two Monte Carlo 
programs are somewhat different from each other. In the string fragmentation picture, given 
an initial q 0 , it is assumed that a new q^qi pair may be created according to some pre-assigned 
probabilities such that a meson q 0 qi is formed and a q\ is left behind. This q\ may at a later 
stage pair off with a g 2 , and so on. The relative probabilities used in PYTHIA to create 
uu, dd, ss, cc pairs are set at 1 : 1 : : y c , where % = 0.3, y c = 10 -11 are default values. 

This choice is based on a quantum mechanical effect and is only approximate. For example, 
in a paper by Lusignoli et al. [|I4|, in order to reproduce the B c meson production using 
HERWIG to within a reasonable range, they had to increase the probabilities of creating cc 
pairs in the cluster splitting. Thus, this ad hoc feature of string fragmentation in these Monte 
Carlo programs makes it less predictive than our perturbative QCD fragmentation functions, 
because in our perturbative QCD fragmentation model the probability for b —> B Sl B* is 
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determined reliably once mb, m s , as well as the nonperturbative parameter of the bound- 
state, e.g., the decay constant /%,, are given. The string fragmentation picture is rather 
phenomenological since the parameter can be adjusted freely to fit the experimental 
measurements. Another weakness of such a string fragmentation model is that it must 
assume the B s /B* ratio of 1/3, from naive spin counting, while our perturbative QCD 
fragmentation model can reliably predict this ratio once m s is given. In fact, the ratio is 
slightly larger (20%) than 1/3, as expected for a non-zero m s , in the HQET JT5[] in which 
m s /mb is the leading heavy quark-spin symmetry breaking effect. 

The comparison with forthcoming data from the Fermilab Tevatron should be instructive. 


This work was supported by the U. S. Department of Energy, Division of High Energy 
Physics, under Grant DE-FG02-91-ER40684 and DE-FG03-93ER40757. 
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TABLES 


TABLE I. Total B s and B* cross sections in nb versus p“ in for pt > 5 GeV and \y\ < 1 at 
\/s = 1.8 TeV. Variation with the scale p are shown. 


^min 

Pt 

v(B s ) 

<B* S ) 

a(B s 




Mi?/ 2 

PR 

2 Mi? 

Mi?/ 2 

PR 

2 Mi? 

Mi?/ 2 

PR 

2 Mi? 

5 

300 

470 

650 

700 

1100 

1500 

1000 

1600 

2200 

10 

45 

62 

65 

110 

145 

150 

150 

210 

220 

15 

14 

16 

14 

32 

37 

33 

46 

52 

47 

20 

5.2 

5.1 

4.4 

12 

12 

10 

17.5 

17 

15 
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FIGURE CAPTIONS 


1. The transverse momentum distributions of B s and B* mesons for p R > 5 GeV, \y\ < 1, 
and n = hr/2,h r , and 2p R at y/s = 1.8 TeV. 

2. The B s and B* total cross sections for p T > p“ in with p T > 5 GeV, \y\ < 1, and 
// = pr/2,pr, and 2y R at yfs = 1.8 TeV. 
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